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A consistent treatment of the ground state correlations be- 
yond the random phase approximation including their influence 
on the pairing and phonon-phonon coupling in nuclei is presented. 
A new general system of nonlinear equations for the quasiparticle 
phonon model (QPM) is derived. It is shown that this system 
contains as a particular case all equations derived for the QPM 
early. New additional Pauli principle corrections resulting in the 
anharmonical shifts of energies of the two-phonon configurations 
are found. A correspondence between the generalized QPM equa- 
tions and the nuclear field theory is discussed. 



1 Introduction 



Many properties of the nuclear vibrational states can be described within 
the Random Phase Approximation (RPA), which enables one to treat some 
correlations in the ground state [1-5]. The low- lying nuclear vibrational 
states investigated with new generation of detectors [H and the double giant 
dipole resonances observed in relativistic heavy ion collision []7|, provide 
an excellent test for the studies of deviations from the harmonic picture 
for multi-phonon excitations in real physical systems. It is well known 
that due to the anharmonicity of vibrations there is a coupling between 
one-phonon and more complex states [§, ^. Usually such a couphng was 
considered for the RPA phonons only [0] . 

From another point of view the RPA violates the Pauli principle and 
many attempts have been done to improve it [9-29] . Extended RPA equa- 
tions include additional corrections for the ground state correlations (GSC) 
n [g, In the papers |2||, ^ phonons of the extended RPA 

containing corrections for the GSC have been used as a basis on which 
the quasiparticle-phonon model (QPM) equations [§, ^ are general- 
ized. The influence of the GSC on properties of nuclear vibrational states 
constructed by one- and two-phonon configurations was studied in [^ . 
Besides the GSC, the Pauli principle corrections arising in the two-phonon 
terms due to the fermion structure of the phonon operators were taken into 
account and the particle-particle channel was included too [^, [29| 



However, these extended QPM equations have been derived for the case 
of the so called quasidiagonal approximation for the Pauli principle correc- 
tions. In the present studies we don't use this approximation. The second 
improvement seems more important. Namely, it is found that there are 
additional Pauli principle corrections resulting in the anharmonicity shifts 
of the two-phonon configuration energies. A correspondence between the 
anharmonicity shifts of the QPM and Nuclear Field Theory diagrams is 
discussed. 



2 Hamiltonian of the model 

We employ the QPM hamiltonian []^, |^ including an average nu- 
clear field described as the Woods-Saxon potential, pairing interactions, 
the isoscalar and isovector particle-hole (p-h) and particle-particle (p-p) 
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residual forces in a separable form with the Bohr-Mottelson radial depen- 
dence 



^ = E f (Ej - Xr)aUa,^ - : (r)Po {r] 

''' \jm ^ 



[(4'"" + ^M''"') : Mlf (r)M<;>(c.r) :] , (1) 

We sum over the proton(]3) and neutron(n) indexes and the notation {r = 
is used and a change r ^ — r means a change p ^ n; a is the 
channel index a = {ph^pp}. The single-particle states are specified by 
the quantum numbers (jm); Ej are the single-particle energies; A,- is the 
chemical potential; G'^^ and k,^^'^ are the strengths in the p-p and in the 
p-h channel, respectively. The monopole pair creation and the multipole 
operators entering the normal products in ([I|) are defined as follows: 



jm 



jj mm 



Af<^>+ (r) = Y. ; {jmj'm | A^,>/Wa+.«t„ ,, 

V I jj mm' 

where /jy^ are the single particle radial matrix elements of residual forces. 
3 Extended RPA 

In what follows we work in the quasiparticle representation, defined by the 
canonical Bogoliubov transformation: 

'^jm ~ '^j^jm "I" {~^y Vjaj^rri- (2) 

The Hamiltonian (|I]) can be represented in terms of bifermion quasiparticle 
operators (and their conjugate ones) [g, |33, |3T|] : 



V A— /i 



B{33 ; A/i) = (-1)-^ ^™ {J^J ^ I V)«/m«j'-m'^ 



^^(j/; A/i) = E {jmj'm I A/i)a+,at^, 



mm 



We introduce the phonon creation operators 



QU = o E A-^Ur, A/i) - ^(jj; A - ^)) . (3) 



jj 



where the index A denotes multipolarity and fi is its z-projection in the 
laboratory system. One assumes that the ground state is the ERPA phonon 
vacuum | 0), i.e. Qx^i | 0) =0. We define the excited states for this 
approximation by Qxf^i \ 0). The following relation can be proved using 
the exact commutators of the fermion operators: 

jj 

where q-' = qj + and qj is the quasiparticle distribution in the ground 
state: qj = Em(0 | a^najm \ 0) {2j + 1)^"^ . The quasiparticle energies 
the energies of the ERPA excited states {co^J, the chemical potentials (A^), 
the coefficients {u,v,ilj,ip) of the Bogoliubov transformations (@), and 
the quasiparticle distributions in the ground state {qj) are determined from 
the following non-linear system of equations obtained by using the general 
equation of motion approach |Tlj]: 

(0| {6aj^, [H, J } I 0) = eM {Sa^^, | 0), (4) 

(0| [6Qx,„ [H, Q+J] I 0) = ujjOl [6Qx,^, QjJ | 0), (5) 
with the constrained conditions: 

{^jm, Oifm'} = ^jm,fm'j (6) 

{0\E^4m^jJO) = N(^\ (7) 

jm 

{0\lQx,.,QtJ\0) =Sii', (8) 

where A^^'^^ is the number of particles. The variations 6ajm and 6Qx^i have 
the following form: 

6ajm = Sujajm - (-l)^""'(5i;jajl^, 
3 
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The system of equations (@)-(0) can be derived using various approaches, 
e.g. [|r?j. For Qj = the system of equations reduces to the usual 

BCS and RPA equations with the p-h and p-p channels |5D| . 
The equation (|^) have been obtained under the assumptions: 



(0| {cyj[m[j [H, ajm]} afm'aj.m, | 0) « 0, (10) 
(0| [H, a+ J aj^m, I 0) « £jqjSjm,j,m,, (11) 

(0| {[H, a+J , ttjimi} aj,^>af_^jn[ \ 0) « £jg/^jmjimi^j'm'j>i- (12) 

If we put Qj = then these approximations correspond to the procedure of 
the linearization of equations. One can derive the equation (g) using vari- 
ous methods [jT2|, [TBj . Using the completeness and orthogonality conditions 
for the phonon operators one can express bifermion operators A'^{jj ; A/i) 
and A(jj'] A/i) as 

A+ijf; A/i) = (1 - ) E ii^jfQU + ^^'Qx-,i) (13) 

i 

It is necessary to point out that the solutions of the system of equations 
0)"(@) obey the following equalities: 

(0| [Qx,.i , [H , Qt,,]] |0) ^^^^ = i^xAi', (14) 
The proof of these statements is analogous to that in the usual RPA case 



As it was shown in Refs.|T7|, |2^, the ERPA calculations give a better 
agreement with experimental data for the characteristics of the low-lying 
states than the RPA ones. 



4 Generalized QPM equations 

The GSC affect not only the RPA, but they also should change 
the quasiparticle-phonon couphng (see Refs.[^, |2D|). To take into ac- 
count such effects we follow the basic ideas of the QPM. Hereafter we 
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generalize the extended QPM equations [j2^, |2^. As it was shown in 



the pairing vibrations give a neghgible contribution to qj. On the 



other hand the two-phonon configurations including the pairing vibration 
phonons have an energy essentially higher than the configurations con- 
structed from usual vibration phonons. That is why we do not take into 
account the coupling with the pairing vibrations (A = 0) in what follows. 

The initial Hamiltonian (|I]) can be rewritten in terms of quasiparticle 
and phonon operators in the following form: 

H = ho + K + hQQ + hQB, (15) 



jm 



hQQ = hi + h2 + hs, 



Xllii' T 

h =-- T 

^ X,^^'r'^Jyr''y- 



yX^yX^' 



Xfii^ Xfii ' 



Xi^Ai 



h,= 



1 ^ Xf (r)+Xf (r) 
8 xkr ^JWW 



QXfiii Q 



Xjii 



hQB 



1 



Q^^L^r> {r))B{jj;X-fi) + h.c., 



The coefficients of the hamiltonian (|T5D are given in Appendix A. The term 
hQB is responsible for the mixing of the configurations. While constructing 
the hamiltonian ( [[51) we have neglected the terms ~: B(jj'; A/i) B(jij[] A/i) : 
which do not lead to coherent effects and the energy corrections due to these 
terms are small in spherical nuclei [j^, ^ . 

The commutation relations Qi,Q2 ^Qt , Qi^Qt ^Qt 

Qsj Qij Qt Q^re calculated by using the transformation (|T3|). 

=E^i(41|23)Q4 + i^(41|23)g+ (16) 

4 
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Q+ Q+] , Q+] = E i^3(41|23)Q4 + K2{41\23)Q+ 

4 

Q-i, [Ql.Qoil = E A'5(41|23)04 + A'4(41|23)Q4+ 



(17) 



(18) 



These coefficients are of fourtfi-order in pfionon amplitudes. In wfiat 
follows one needs the coefficients Ki, K4 only. They are given in 
Appendix B. We obtain the anharmonic corrections taking into account 
these commutation relations. 

Hereinafter the ground state is approximated by the ERPA phonon 
vacuum. Using the ERPA phonons as a basis the wave functions of the 
excited states of even-even nuclei can be written as: 



where 



(19) 



i AiiiA2i2 

Q \lfililQ\2^^2i2) \n — ^ ^>^U^l>^2^^2Q\WlilQ\2m2- 



The wave functions (|TD|) have the normalization condition 

>2 



AiiiA2i2 Xiliil X2'^2' 

. .A1+A2+J 

A2«2 



(20) 



(^A^z^,A2i2<^A'iz'i,Ain + ^A^i^,Ain<^A'ii'i,A2i2 (-1)^^+^'+'^ + 



K-^{X2'l2,Kh I Ai?i,A2^2)) = 1- 

We use the equation of motion method to diagonalize the hamiltonian ([T5|) 



(o| [sn 



0) = E.(0|[(5Qa,.,^^V] |0), 



where the variation 5^\^y has the following form 
8^y^,, = Y.5R,{\v)Qx,^+ E 5Ptl'M^){Q^2,2^2Qx^ 



XiiiX2i2 



'X2H2i2^MfJ.iii)Xfi 



1 



,Ai+A2+A 



The variational principle yields a set of linear equations for the unknown 
wave function coefficients Ri{Ju) and P^\*^(Jz/) 



I 
X 



Viv) 



'2V 



with the additional condition (|2tl|). The matrix lA^ is a transpose of the 
matrix lA. The number of linear equations ([21| ) equals to the number of 
one- and two-phonon configurations included in the wave function 
The notations introduced above are 



Wif (A) = (o|[gA,.,[^,Qj,J]|o). 



(22) 



X' i' y i' 

2 Aiii,A2j2 



(A) = (-1) 



X 



(23) 



(0| [{Q\'^t.'AQ>^w)x^,^ (^AiMi^i^ 



- A2M2j2yA/x 



|0>. 



|0), 



(24) 



^^l'lj'^2*2 t\\ 

Aiii,A2i2 ^ ' 



A1+A2+A 



+ 



h'^i'^Mi2^KAMii + ^A^z;,,Aizi<^A'iii,A2i2 
A'\A'/2,AX I AiZi,A2Z2)). 

Now we calculate the matrix elements (^), (|2^) and as a result we get 



X 



(25) 



Z^"' (r) (-1)^+^1+"^' [K(\\ii'^, Xm I A^, - Ki'iXii'^, Xi \ Xm, Xi')) + 
Xf (r) + Xf * (r)) A'^'(Ai^'2, | Xi\ Xm) 



>^a::;;aS (a) = {s^',,, mf, (ao + 5,^,, {x,)) + 

<5ai,A2^a^,a, (-1)^^+^^+^ {5,^,^ (A2) + 5,^,, (Ai)) + 

Aa;^(A2i2, Ai^i | AiH,A2i2)- 

have been evaluated by keeping all terms 



The matrix elements 
containing the first power of coefficients Ki and products KK^ only. We 
took into account also the terms containing the products KKi that include 



the ijj'^^p terms. The other terms including the products KKi are smaller 
and they are neglected. The terms with products KiKi were not taken into 
account for these matrix elements because they have no i/j^^ip^cp terms. It 
should be pointed out that the matrix element ( |25| ) differs from the phonon 
energy in contrast with the ERPA case (see Eg . ) . This difference is 
due to the approximations (PD|) , (\n\) , (P^) . The corrections to the matrix 
element include phonon amplitudes of the next orders. Only the terms 
/i2 and /i3 of the hamiltonian ([TBD contribute to these corrections. Then, 
these corrections become zero cp = (see the explicit form for Xi, Ki, 
K2). It should be mentioned that we do not use so called quasidiagonal 
approximation for the coefficients K and Ki. 

The values Acj^ are anharmonic shifts of energies of the two-phonon 
configurations due to the Pauli principle corrections, where AcJ^ = /SjJj\ + 
AcJs + AcJ^, 

/^u\{\'2i'2, Kh I AiH,A2^2) = (26) 

(cJAiii + Ux^i^) K^{>^2H^KH I Alii, A2i2)- 




K (A2^2, Ai^ I Ai^i, A2^2) + 



— A (Ai^i, I Ain, A2^2) (-1) + 




it'''' (r) + X^'^' (r) 



X 



K^iX-iiz.X^i I Ain, A2^2)i^^(A2^2, A^ii | 



/\ul{M2H^KA I AiH,A2^2) = 



1 _ (Z^^'^' (r) + 



X 



(27) 



A'i^(A2^2, Ai^i I X2i,X'ii'i 



(r) + (r) 



X 



A'2(A2i2' A'lZ I AiH, A2^2) 



Z>^2^2^ (r) + Z 



X 
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(-l)^"'''^+'^i^2(A;z;,A'2Z I AiZi,A2i2)- 

.BA-^yPy^^ (A2Z2,M|A.Z,A3Z3) 

i^2\A3Z3, A4Z' I AiZi, A2Z2) " E (- X 

J1J2 



i^"^'(A3i3, AlZi I A4Z, Aii'i) X 



X 



/ 


' Jl 


A3 


A.| 




A4 


J2 


A2 


v 


. -^1 


A'2 


A 



Ki\>^A. Mi2 I A4i', A3i3) + (-1)^3+Ai+J.+J.+A ^ 

i^^^(A3Z3,AlZl I A4Z,A'/2)A^/^(A3Z3,A2Z2 | , KA))) + 



Jl A3 Ai 
A4 J2 A2 
A2 A'l A 



^^A.+A3+A.+A (2A3 + 1) ^^^--^ (r) 



(2Ai + 1) 



X 



i^2^^(Al23, A4i I Alii, A40^\A'/2, A'li'i I A2Z2, Aii3) + (-1) 



a/ 2Ai + r 

,2A2 + 1. 



i^2^(A2Z3,A4i I A2Z2,A40i^'(A'2z'2,A;z; I AiZi,A2Z3)))) 



X 



Aa;^(A'2i2, A'li'i | Aiii, A2i2) = -7 E 



i^4^^(A'2i'2,A2i| AiZi,A;i;)(-l) 



3 V 



vA+A3+A2+Ai+l 



r) + X"^"^ (r) 



X 



(2A3 + 1) 



Ai A2 A 
A2 A'l A3 



1 XY'' (r) + XY'' (r) 



'2A3 + 1 



(28) 



+ 



7V^ M^^HAli4, A4Z I A4i', AiZi)A'\A'2i'2, A'li'i | Aii4, A2Z2) + 
(^^^) ^^'(^2^4, A4i I A4i', A2i2)i^\A'2z'2, A'li'i | Ai^i, A2i4))) • 

One can prove that AcJ^ ;> Aa;2 and AoJ^ AoJg in the case ip ^ (f. 
Moreover, the shifts Aa;2 and AcJ^ become zero if (/? = 0. It is necessary 
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to emphasize that only the terms /i2 and contribute to the shifts AlJ2 
and Aa;3. The shifts Aa;2 and Aa;3 have been neglected in our previous 
papers. The matrix elements coupling one- and two-phonon configurations 
(^) are 



T 



where 



C^iA2{Mr) = (-l)*'+^=+*^\/(2Ai + l)(2A2 + l) £ '(l-fe)^ 



J1J2J3 



^2tS (A^^ r) = (-l)^^+^^+^4^V(2Ai + 1)(2A2 + 1) E ^(l-^,.,3)x 



72 



J1J2J3 

:^^2*2 \ ji j3 j2 



^^A;;;(A^,r) = E (^2 03 (A^,r)A^\AYA3Z3 | Ai^i, A2^2)+ 

A'i'A3l3 

'A/ 



2A' + 1 

/73^;X (A^, r)A'^(A2^2AY|A3Z3AiZi, ) + E (2/ + 1) ^ ^X+r"" ^^^^ 
} {k'\X'i'Ai I Aaza, A2Z2) ^i^l^^// (A3Z3,r) (-i)^+^'+^^+^3 + 
A^/'(AYA2^2 I A3^3, Az) U,^i^ (A3^3,r) (-1)^+^'+-^' 



I (i^^'(AY, Az I A3.3, Ai^i) (A3^3,r)+ 
A^/'(AYAiZi I X,z,,Xi) U,^i (A3^3,r) (-l)^^+^-3+^' 
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72 



J1J2J3 




J1J2J3 




V, ^2^2 Alii , £Ai I \ My \2i2\ 



The terms including the coefficients Ki in were neglected in the papers 
21, iZI, |2g. 



Solving the system of equations ( |7T|) we obtain the energies {E^y) of ex- 



cited states and the coefficients Ri{Ju) and P^^l^{Ju) of the wave function 



These equations are more general than ones derived in [j2^ |27], 



They have the same form as the basic QPM equations []5|, The GSC 
change phonon energies u\i^ the anharmonic shifts of two-phonon configu- 
rations and the matrix elements coupling one- and two-phonon configura- 
tions. In the case when i/j ':$> (f the system of equations ( |2T|) can be reduced 
to the system of equations of the extended QPM |27|, If we put 
Qj = 0, Xi = and i^'i = we get the usual QPM equations with taking 
into account the Pauli principle corrections []^, |^ In the case when 
Qj = 0, K = and Ki = we have equations describing coupling of one- 
and two- RPA phonons without taking into account the Pauli principle . 

5 Anharmonic shifts of two-phonon configurations 

Now we discuss the correspondence between the QPM equations presented 
above and the diagrams of the Nuclear Field Theory (NET) For 
the sake of simplicity we compare these approaches for the hamiltonian ([I|) 
including the average nuclear field and the isoscalar particle-hole residual 
forces only. Besides that, we put Qj = in the system of equations (|21|), 
i.e. conventional RPA phonons are used as the QPM basis in this case. 

The NFT is a formulation of many-body perturbation theory with vi- 
brational modes summed to all orders in RPA. Its building blocks are RPA 
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phonons and the single particle degrees of freedom which are described in 
the Hartree-Fock approximation. The coupling between them is treated di- 
agramatically in the perturbation theory. Diagrams illustrating first order 
coupling between the surface vibrations and the fermion fields are shown 
in Fig. 1. The wavy lines are phonon propagations, while the particles and 
the holes are depicted by the arrowed lines. The lowest-order anharmonic 
terms of NFT contributing to the energy of two-phonon states are repre- 
sented by fourth-order diagrams shown in Fig. 2 [3,35-41]. These graphs 
are called butterfly-type (A,B), trapezoid- type (C,D), and diamond-type 
(E,F) diagrams. For each diagram shown in Fig. 2 there are 5 other dia- 
grams which are obtained by changing the direction of the phonon lines. 

One can rewrite the system of equations (^) in the space of two-phonon 
states. The diagonal approximation for the matrix element Wi* (A) 
( Wi • (A) = 6i^i' Wi ■ (A)) enables one to find the coefficient Ri{Jv) from the 
first equation of the system (^. Substituting it into the second equation 
of this system one can get the following secular equation to find energies 

det II (^,.,, Wif; (Ai) + 5,^,, Wii (A2)) ^Ai,A,^A'„A.+ 
h',,xMM (-1)'^+'^+' (<5,^,, (A2) + 5,^,, (Ai)) + 

AcJ^(A2^2, I Ain, A2^2) - At/^(A2^2, Kh I Ain, A2^2) - EJ\ = 0, 
where 

At/ (A2^2, Ai^i I Am, A2^2) = E ..J ^(.. ' • (30) 

It is seen that the values At/ and AcJ are anharmonic shifts of two- 
phonon states. The term Kqb does not contribute to the energy shifts of 
AcJ, while it results in the energy shifts of At/. Furthermore, the term Yiqe 
contributes to the scattering vertices (see Fig.l (A), (B)) only. A proof of 
this statement is similar to one of [j3Bj. For example we consider terms /S.K 
(see (g6|)) and At/2 (see (0)): 

AA'^(A4i4, As^a I Ain,A2i2) = 

-7 EE . , ^ (A4^4,A3^5 I AiZi,A2^2), 



«5 ''•'1^5 
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The anharmonic shifts of energies AK {X4i4, X^i^ \ Ai^i, A2^2) and 
AU2{X3i3, X^i^ I A2^2, Aih) can be illustrated by the QPM diagrams [§, ^ 
shown in Fig. 3 (A) and (B), respectively. If the intermediate phonon 
line (^5) in these shifts are changed by the two-quasiparticle state lines 
(i.e. il) = 1 and ip = 0), then these shifts can be represented by the 
fourth-order diagrams of the NFT (Fig. 2). In this case the term AK 
corresponds to the butterfly-type diagrams, while the term AU2 can be 
presented as the diagrams of the trapezoid-type and the diamond-type. 
Using the same method one can prove that the shifts Au and AU (see 
( P5H28| , [3D|) ) correspond to diagrams shown in Fig. 2 (A), (B) and Fig. 2 
(C), (D), (E), (F), respectively. The shifts Aa;2 and Azj^ correspond to the 
butterfly-type diagrams. 

6 Conclusion 

A consistent treatment of the ground state correlations beyond the RPA 
including their influence on the pairing and the phonon-phonon coupling 
in nuclei is presented. A new general system of nonlinear equations for 
the quasiparticle phonon model is derived. It is demonstrated that this 
system contains as a particular case all equations derived for the QPM 
early. The new additional Pauli principle corrections resulting in the an- 
harmonic shifts of energies of the two-phonon configurations are found. It 
is shown that the anharmonic shifts due to the Pauli principle corrections 
correspond with butterfly-type diagrams of the nuclear field theory, while 
the anharmonic shifts due to the matrix elements coupling one- and two- 
phonon configurations give rise to the trapezoid-type and diamond-type 
diagrams of the NFT. 
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Appendix A 

The coefficients of ttie tiamiltonian ([T^) are given by tfie following expres- 
sions: 



X^'' (t) 



^D^' (r) + D'l' (r) (r) + (r) (r 



^'^^ (r) = ^ ^ (Do'^ (r) (1 - t^' (r)) + F,'^ (r) + (D^ (r) + D>^^ (r 
z^^' (r) - 4^' (r)) + Fi^^ (r) (r) - F^^^ (r) (r)) + 



Xi 



X 



^ ^ (Do'^' (r) (1 + to'^ (r)) - Fo^^' (r) + (d^^' (r) - D^' (r) ) x 
z'^ (r) + 4^ (r)) + (r) F^^^' (r) - (r) F,^^' (r 



Xf ^ (r) 



M ^ (^0^ ir) (2 - to'' ir)) " i^o'^ (r) + {D>^' (r) - Z^? (r)) x 



z^^ (r) - (r) ) + tf (r) F^^ (r) - tf (r) F,'^ (r 



/-Ai' zT'Ai 



l'^^^) (r) = ^^7,^ + (z^^ (r) ± (r 



(-) , ^^^^ 



u^7) - 

3] JJ 



2(2A + 1)^ 



D^^ (r) + /s:S''^')) + Do'^ (-r) f^^'^^^^ - ^(^'^'^^V 



Dil' (r) (/^E)"'^^^ + «:S''^^^) + (-r) (4"'^^^ - 4''''''^ 



n = {+,-} 
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"'"''^ D^^ (r) (4^'^') + A^l''^')) + D^^ (-r) (4^'^^^ - a^I^'^'^ 



ni = {0,l,2} 



. .1 

33 



33 



33 



33 

(±) I (±) I 

33 ■' ■' ■' ■' 33 ■'3 ■'3 

Appendix B 

The coefficients K have the following form: 

K {X2fi2i2^ X[fl[i[ I Ai/iiZi, A2/i2^2), 



K^iX'^i'^A'^' I Az,A2Z2) = V(2A + 1) (2A2 + 1) (2A' + 1) (2A'2 + 1) x 





( 




h 


aU 




< 


k 


h 


A' ^ 


jlj2j3j4 




A 


A2 


J 



15 



^•^^1^2 21,31 ^J3J4 ^Jlj4 Jij2 "t^JSjl Jaji'f Jij4 J ^ 

(_1)J1+J2+J3+J4+A+A2+A'+A^+J I -^2 I J K X 
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Figure 1: Diagrams illustrating the first order coupling between the surface vibrations 
(wavy lines) and the fermion fields (arrowed lines). 




(A) (B) 



Figure 3: The QPM diagrams (A) and (B) illustrating the anharmonic shifts 
AK^{X4^i4, X^i^ I Al^l,A2^2) and A[/2 (-^3^3) ^^4^4 I ^2'i'2,Xiii) of energies, respectively. 
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